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The Ni ion in LaNi02 has the same formal ionic configuration id? as does Cu in isostructural 
CaCu02, but it is reported to be nonmagnetic and probably metaUic whereas CaCu02 is a magnetic 
insulator. From ab initio calculations we trace its individualistic behavior to (1) reduced 3d — 2p 
mixing due to an increase of the separation of site energies (e^ — Sp) of at least 2 eV, and (2) 
important Ni 3d(3z^ ~ mixing with La 5d(3z^ ~ t"^) states that leads to Fermi surface pockets 
of La 5d character that hole-dope the Ni 3d band. Correlation effects do not appear to be large in 
LaNi02. However, ad hoc increase of the intraatomic repulsion on the Ni site (using the LDA+U 
method) is found to lead to a novel correlated state: (i) the transition metal d{x^—y^) and d(3z^ — r^) 
states undergo consecutive Mott transitions, (ii) their moments are antialigned leading (ideally) to 
a "singlet" ion in which there are two polarized orbitals, and (iii) mixing of the upper Hubbard 
3d{3z^ — r^) band with the La 5d{xy) states leaves considerable transition metal 3d character in a 
band pinned to the Fermi level. The magnetic configuration is more indicative of a Ni^"*" ion in this 
limit, although the actual charge changes little with U. 

PACS numbers: 71.20.-b, 71.20.Be, 71.20.Eh, 71.27.+a 



I. INTRODUCTION 

The perovskite oxide LaNiOs, purportedly an ex- 
ample of a correlated metallic Ni'^^ system, has been 
investigated over some decades by a few groups 0, 
0, 13 for possible exotic behavior. The oxygen-poor 
lanthanum nickelate LaNiO^; has also attracted at- 
tention, because of characteristic changes of its elec- 
tronic and magnetic properties as the oxygens are 
removed. It is metallic at 2.75 < x < 3, but semi- 
conducting for 2.50 < X < 2.65. 4] For x — 2.6, 
it shows ferromagnetic ordering with 1.7 /is/Ni be- 
low 230 K y and magnetic behavior of the x — 2.7 
material has been interpreted in terms of a model 
of ferromagnetic clusters. jj| At x = 2.5, where for- 
mally the Ni is divalent, a perovskite- type compound 
La2Ni205 forms in which NiOe octahedra lie along 
c axis directed chains and Ni04 square-planar units 
alternate in the a — b plane. This compound shows 
antiferromagnetic ordering of the NiOe units along 
the c axis but no magnetic ordering of the Ni04 
units. 

Since LaNi02 with formally monovalent Ni ions 
was synthesized by Crespin et a/.0, it has at- 
tracted interest0,[i3,[i3 because it is isostructural 
to CaCu02,^3 the parent "infinite layer" material 
of high Tc superconductors, and like CaCu02 has 
a formal cfi ion amongst closed ionic shells. How- 
ever, it is difficult to synthesize and was not revis- 
ited experimentally until recently by Hayward et al. 
who produced it as the major phase by oxygen dein- 
tercalation from LaNiOs. 13] Their materials consist 
of two phases, the majority being the infinite-layer 
(Ni02-La-Ni02) structure and the minority being a 
disordered derivative phase. Magnetization and neu- 



tron powder diffraction reveal no long-range mag- 
netic order in their materials. Its paramagnetic sus- 
ceptibility has been fit by a Curie- Weiss form in the 
150 < T/K < 300 range with S=i and Weiss con- 
stant 9 — -257 K, but its low T behavior varies 
strongly from this form. More recently, this same 
group has produced the isostructural and isovalent 
nickelate NdNi02.[l3 

One of the most striking features of LaNi02 is that 
it potentially provides a structurally simple example 
of a monovalent open shell transition metal (jP ion. 
Except for the divalent Cu^+ ion, the d^ configura- 
tion is practically nonexistent in ionic solids. In par- 
ticular, the formal similarity of Ni^+ and Cu^+ sug- 
gests that Ni^+ compounds might provide a "plat- 
form" for additional high temperature superconduc- 
tors. It is these and related questions that we ad- 
dress here. 

In this paper we present results of theoretical 
studies of the electronic and magnetic structures 
of LaNi02, and compare with the case of CaCu02 
(or isovalent Cai_a;Sra;Cu02) which is well charac- 
terized. A central question in transition metal ox- 
ides is the role of correlation effects, which are cer- 
tainly not known a priori in LaNi02. We look at 
results both from the local density approximation 
(LDA) and its magnetic generalization, and then ap- 
ply also the LDA+U correlated electron band theory 
that accounts in a self-consistent mean-field way for 
Hubbard-like intraatomic repulsion characterized by 
the strength U. Our results reveal very different be- 
havior between LaNi02 and CaCu02, in spite of the 
structural and formal (P charge similarities. The dif- 
ferences can be traced to (1) the difference in id site 
energy between Ni and Cu relative to that of Cu, (2) 
the ionic charge difference between Ca^"*" and La'^"'" 




FIG. 1: (Color online) Crystal structure of LaNi02, 
isostructural to CaCu02. Ni ions are in the origin and 
La ions in the center of the unit cell. It has no axial 
oxygens. 



and associated Madelung potential shifts, and (3) 
the participation of cation 5d states in LaNi02. 

We also discuss briefly our discovery of anoma- 
lous behavior in the transition metal 3(P ion as de- 
scribed by LDA-I-U at large U. Although well be- 
yond the physical range of U for LaNi02, we find 
that LDA-I-U produces what might be characterized 
as a "singlet" ion in which the internal configu- 
ration is one d{x'^ — y^) hole with spin up and one 
d{?tz^ — r^) hole with spin down, corresponding to 
an extreme spin-density anisotropy on the transition 
metal ion but (nearly) vanishing net moment. 



II. STRUCTURE AND CALCULATION 

In the samples of LaNi02 synthesized and re- 
ported by Hayward et al., there exist two phases 
with space group PA/mmm (No. 123) but different 
site symmetry. We focus on the majority infinite- 
layer phase, which is isostructural with CaCu02. 1^ 
In the crystal structure shown in Fig. 1, Ni ions 
are at the corners of the square and La ions lie at 
the center of unit cell. The bond length of Ni-0 
is 1.979 A, about 2% more than that of Cu-0 in 
CaCu02 (1.93 A). We used the lattice constants 
a = 3.870931, c = 3.3745A,0 with a {^2 x V2) 
supercell space group lA/mmm (No. 139) for AFM 
calculations. 

The calculations were carried out with the 
full-potential nonorthogonal local-orbital (FPLO) 
mcthodllSl and a regular mesh containing 196 k 
points in the irreducible wedge of the Brillouin 
zone. Valence orbitals for the basis set were 
La 3s3p3d4s4p4d5s5p6s6p5d4/, Ni 3s3p4s4p3d, O 
2s2p3s3p3(i. As frequently done when studying tran- 
sition metal oxides, we have tried both of the popu- 
lar forms of functional [l^fl7l | of LDA-I-U method [Tsll 
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FIG. 2: LDA paramagnetic band structure of LaNi02. 
The Ni 3d{x^ — y^) band crosses the Fermi level (zero 
energy) very much as occurs in cuprates (see Fig. 3). 
The La 4f bands lie on 2.5 - 3.0 eV. The La 5d(3z^ - r^) 
band drops below Ef at F and A. 

with a wide range of on-site Coulomb interaction U 
from 1 to 8 eV, but the intra-atomic exchange inte- 
gral J=l eV was left unchanged. For CaCu02, we 
used the same conditions as the previous calculation 
done by Eschrig et al. using FPL0.[l9j 



III. RESULTS 
A. LDA Results 

We present first the LDA results. The param- 
agnetic (PM) band structure with its energy scale 
relative to Fermi energy Ep is given in Fig. 2. A 
complex of La 4/ bands is located at -1-2.5 eV with 
bandwidth less than 1 eV. The O 2p bands extend 
from about -8 eV to -3.2 eV. The Ni 3d bands are 
distributed from -3 eV to 2 eV, with the localized 
t2g complex near -1.5 eV, while the broad La 5d 
states range from -0.2 eV to 8 eV. Unlike in PM 
CaCu02, there are two bands crossing Ep. One is 
like the canonical d{x'^ — V^) derived band in the 
cuprates, rather broad due to the strong dpa an- 
tibonding interaction with oxygen Px,Py states and 
enclosing holes centered at the M point. The other 
band, lying at -0.2 eV at F and also having its max- 
imum at the M=(-, -,0) point, is a mixture of La 
5d{3z^ — r^) states and some Ni 3d{3z^ — r^) char- 
acter. Already this band indicates importance of Ni 
3(i - La bd band mixing. 

Using a simple one-band tight binding model 

£k= £o — ^ tfi e^'^'^, 

R 



2 




-2.0 I I I I I I I 

rxM rzRA z 

FIG. 3: "Fatband" representation of Ni M(x'^ — y^) in 
LDA. This band appears at first very two-dimensional, 
but is not because (1) the saddle point at X(0,7r/a,0) is 
not midway between the F and M(7r/a, 7r/a, 0) energies, 
and (2) kz dispersion between the X and R(0, vr/a, vr/c). 



TABLE I: Tight binding parameters (in meV) for Ni 
3d(x^ - ) of LaNi02 and Cu 3d(K^ - y ^ ) of CaCu02 . eo 
is the site energy and t's are hopping integrals. Ratio (in 
%) is hopping integrals for LaNi02 to those for CaCu02. 



parameters LaNi02 CaCu02 \Ratio\ (%) 



eo 


93 


-200 




t(lOO) 


381 


534 


71 


t(llQ) 


-81 


-84 


96 


t(OOl) 


58 


83 


70 


t(lOl) 





-2 







-14 


-19 


74 



the Ni 3d{x^—y^) band shown in Fig. 3 can be repro- 
duced with a few hopping amplitudes, but requiring 
more than might have been anticipated. The site en- 
ergy is So = 93 meV, slightly above the Fermi level, 
and the hopping integrals (in meV) are t(lOO) = 381, 
^(110) = -81, ^(001) = 58 and i(lll) = -14. There 
is no hopping along the (101) direction. As antici- 
pated from the cuprates, the largest hopping is via 
i(lOO). However, to correctly describe the kz disper- 
sion from X-R {i.e. along t: / a,Q^kz) together with 
the lack of dispersion from F - Z (0, 0, /cz) and also M- 
A (tt/o, 7r/a, fcz), the third neighbor hopping terms 
t(lll) must be included. 

The comparison of the single band tight binding 
parameters with those of CaCu02 is given in Table 
I. It should be noted that the state in mind is an 

— symmetry state that is orthogonal to those 
on neighboring Ni/Cu ions, i.e. an — sym- 
metry Wannier orbital. In Ni, the on-site energy 
is 0.3 eV above what it is in CaCu02, lying above 
Yip rather than below. This difference is partially 
due to the different Madelung potential in the two 




FIG. 4: (Color online) Paramagnetic Fermi surface in 
the local density approximation. In the center (not visi- 
ble), i.e. F, there is a sphere (a radius 0.25(7r/a)) having 
d{2,z^ — r^) character of Ni and La. The cylinder with 
radius 0.8(7r/a) contains Ni d{x^ — y'^) holes, whereas an- 
other sphere (a radius 0.4(7r/a)) at each corner contains 
Ni d{zx) electrons. 



differently-charged compounds, but it also reflects 
some intrinsic hole-doping in the nickelate that leads 
to a lower Fermi level. The largest hopping am- 
plitude (the conventional t) is 71% of its value in 
the cuprate, while the second (i') is essentially the 
same. The t(OOl) = tz is also 70% of its value in the 
cuprate, while the other amplitudes are the almost 
unchanged. 

The LDA Fermi surfaces are shown in Fig. 4. As 
for the cuprates, the Fermi surface is dominated by 
the M-centered hole barrel. In this system neighbor- 
ing barrels touch at R—{Tr/a, 0, tt/c) because the sad- 
dle point at R happens to lie at E^?. The Fermi sur- 
faces also include two spheres containing electrons. 
The sphere at F, with mixed Ni and La d{3z^ — r^) 
character, contains about 0.02 electrons. The A- 
centered sphere is mainly Ni d(zx) in character and 
contains approximately 0.07 electrons per Ni. The 
barrel, whose radius of 0.8 tt/o in the (1,1, fc^) di- 
rection is almost independent of kz but which varies 
along (1,0,^:^), possesses about 1.1 holes, account- 
ing for the total of the 1.0 hole that is required by 
Luttinger's theorem and also fits the formal Ni^+ 
valence (which, being a metal and also mixing with 
La as well as with O states, is not very relevant). 

To investigate magnetic tendencies, attempts to 
find both ferromagnetic (FM) and antiferromagnetic 
(AFM) states were made. A stable %/2 x ^2 AFM 
state was obtained, with spin moment 0.53 /is per 
Ni. This state has lower energy by 6 meV/Ni than 
that of PM state. Just as for the paramagnetic case, 
the AFM state has entangled bands of La 5d, Ni 3d 
and O 2p character near the Fermi energy. In con- 
trast to the unpolarized case (and CaCu02), with 
AFM order the large electron pocket has primarily 
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FIG. 5: LDA antiferromagnetic band structure of 
LaNi02. The Ni 3d bands lie above -3 eV and are disjoint 
from the O 2p bands (not shown) which begin just below 
-3 eV. The antiferromagnetism introduces the gap in the 
Ni dpa band midway between F and M in the range 
- 1 eV. The symmetry points are given such as (0,0,x) 
for F(Z) , (l/2,l/2,x) for X(R) and (l,0,x) for M(A). x 
is zero for the first symbols and 1 for the symbols in 
parentheses. 



La 5d{xy) character and the shghtly occupied elec- 
tron pocket at F has a combination of La 5d{3z^—r^ ) 
and Ni 3d{3z^ — r^) character. Attempts to obtain 
a FM solution always led to a vanishing moment. 



B. Consideration of Correlation with LDA+U 

As noted in the introduction, no magnetic order 
has been observed in LaNi02, either by magneti- 
zation or by neutron scattering. Although the lo- 
cal density approximation often does quite well in 
predicting magnetic moments, for weakly or nearly 
magnetic systems renormalization by spin fluctua- 
tions becomes important pol l2ll |23 | and such effects 
are not included in the local density approximation. 
There is also the question of the strength of cor- 
relation effects due to an intra-atomic repulsion U 
on the Ni site. Analogy to CaCu02 (same formal 
configuration, neighboring ion in the periodic table) , 
which is a strong antiferromagnetic insulator, sug- 
gests that effects due to U might have some impor- 
tance. Here we apply the LDA-I-U "correlated band 
theory" method to assess effects of intra-atomic re- 
pulsion and compare with observed behavior, fn the 
following subsection we compare and contrast with 
CaCuOa. 

Upon increasing U from zero in the antiferromag- 
netically ordered phase, the spin magnetic moment 
of Ni increases from the LDA value of 0.53 to 



a maximum of 0.8 /is at U — 3 eV. Surprisingly, 
for C/ > 4 eV the moment steadily decreases and by 
U — 8 eV it has dropped to 0.2 /is/Ni, which is less 
than half of its LDA value. We emphasize that this 
behavior is unrelated to the observed behavior of 
LaNiOa (which may need little or no additional cor- 
relation beyond LDA). However, this unprecedented 
response of the transition metal ion to the imposi- 
tion of a large U gives new insight into a feature 
of the LDA-I-U method that has not been observed 
previously. 

This "quenching" of the local moment with in- 
creasing U results from behavior of Ni 3d(3z^ — r^) 
states that is analogous to those of the 3d{x'^ ^ y"^), 
but with the direction of spin inverted (then with 
additional complications). As usual for a d^ ion in 
this environment, the majority 3d{x^ — y^) state of 
Ni is fully occupied even at J7 = 2 eV, while the mi- 
nority state is completely unoccupied at J7 = 3 eV, 
where the moment is maximum and the system is 
essentially Ni^+ S=^. One can characterize this sit- 
uation as a Mott insulating 3(i(x^ — y^) orbital, as in 
the undoped cuprates. ki U — i eV, the density of 
states has a quasi one-dimensional van Hove singu- 
larity due to a flat band just below (bordering) the 
Fermi energy as can be seen in the DOS shown in 
Fig. El Upon increasing [/ to 4 eV, rather than re- 
inforcing the S — ^2 configuration of Ni and thereby 
forcing the La and O ions to cope with electron/hole 
doping, the Ni d{3z^ — r^) states begin to polarize. 
The charge on the Ni ion drops somewhat, moving 
it in the Ni^+ Ni^+ direction, with the charge go- 
ing into the La 5c? - O 2p states. Idealizing a bit, 
one might characterize the movement of (unoccu- 
pied) majority character of 3d(?tz^ — r^) well above 
Ei? as a Mott transition of these orbitals, which is 
not only distinct from that of the 3)d{x^ — y^) states, 
but is oppositely directed, leading to an on-site "sin- 
glet" type of cancellation. 

This movement of states with increasing U has 
been emphasized in Fig. for easier visualization. 
The resulting spin density on the transition metal 
ion at [/ = 8 eV is pictured in Fig. There is strong 
polarization in all directions from the core except for 
the position of nodes. The polarization is strongly 
positive (majority) in the lobes of the 'id{x^ — y^) or- 
bital, and just as strongly negative (minority spin) 
in the lobes of the 3d(3z^ — r^) orbital. The net 
moment is (nearly) vanishing, but this results from 
a singlet combination (as nearly as it can be repre- 
sented within classical spin picture) of spin-half up 
in one orbital and spin-half down in another orbital 
that violates Hund's first rule. The magnetization 
density is large throughout the ion, but integrates 
to (nearly) zero. 

This behavior is however more complicated than 
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FIG. 6: (Color online)Change of the Ni 3d(3z^ - r^) and 
3ci(a;^ ~ 2/^) densities of states as on-site Coulomb inter- 
action U increases. One can easily identify a splitting 
("Mott transition") of the 3ci(a;^ — y^) states occurring 
near U =0, and the light (green) lines outline their path 
with increasing U (majority is solid, minority is dashed). 
A distinct Mott transition involving oppositely directed 
moment of the 3ci(32:'^ — r^) states is outlined with the 
dark (purple) lines. This moment is oppositely directed. 
The conceptual picture is also complicated by the split- 
ting even at U=Q which persists in the majority states, 
leaving a band at Ef with strong Ni 3d(32^ ~ f^) char- 
acter as well as the expected upper Hubbard band at 4 
eV. 

a Mott splitting of occupied and unoccupied state, 
as can be seen from the substantial Ni id character 
that remains, even for U — 8 eV, in a band strad- 
dling Ep while the rest of the weight moves to ~4 
eV. In both of these bands there is strong mixing 
with La 5d{xy) states. What happens is that as the 
"upper Hubbard 3d{3z^ — r^) band" rises as U is 
increased, it progressively mixes more strongly with 
the La 5d{xy) states, forming a bonding band and 
an antibonding band. While the antibonding com- 
bination continues to move upward with increasing 
[/, the bonding combination forms a half- filled band 
which remains at Ep. 

Thus we have found that for the Ni^+ ion in 
this environment, increasing U (well beyond what 
is physically plausible for LaNi02) results in S" = i 
Ni"'""'^ being converted into a nominal Ni''"^ ion (the 
actual charge changes little, however) in which the 
two holes are coupled into an intraatomic 
singlet. This behavior involves yet a new kind of 




FIG. 7: Isocontour plot of the spin density of the "sin- 
glet" Ni ion ([/ = 8 eV) when there is an — hole 
with spin up and a 3z^ — hole with spin down. Dark 
and light surfaces denote isocontours of equal magnitude 
but opposite sign. 

correlation between the 3d(3z^ — r^) states and the 
3d{x^ — y^) states, but one which is due to (driven 
by) the local environment. 

This behavior is quite different from the results 
for U =8 eV reported by Anisimov, Bukhvalov and 
Rice[l3 using the Stuttgart TBLMTO-47 code. 
They obtained an AFM insulating solution analo- 
gous to that obtained for CaCu02,|JJ.J with a sin- 
gle hole in the 3d shell occupying the 3d{x^ — y^) 
orbital that antibonds with the neighboring oxy- 
gen 2pa- orbital. The reason for this different result 
is not known, but it is now well established that 
multiple solutions to the LDA-I-U equations often 
exist.lHlll 



IV. COMPARISON WITH CaCu02 AND 
DISCUSSION 

Although Ni+^ is isoelectronic to Cu+^, both the 
observed and the calculated behavior of LaNi02 
are very different from CaCu02. In contrast to 
CaCu02, LaNi02 is (apparently) metallic, with 
no experimental evidence of magnetic ordering for 
LaNi02. The differing electronic and magnetic prop- 
erties mainly arise from two factors. First, the Ca 3d 
bands lying in the range of 4 eV and 9 eV are very 
differently distributed from the broader and lower 
La 5d bands in the range of -0.2 eV and 8 eV. Sec- 
ondly, in CaCu02, O 2p states extend to Fermi level 
and overlap strongly with Cu 3d states, and the dif- 
ference of the two centers is less than 1 eV, as can 
be seen in Fig. |S| Thus, there is a strong 2p — 3d hy- 
bridization that has been heavily discussed in high 
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FIG. 8: (Color online) Comparison of LDA projected 
paramagnetic DOS LaNi02 (upper panel) and CaCu02 
(lower panel). Note the separation of the Ni 3d states 
from the O 2p states in the upper panel, which does not 
occur for the more strongly hybridized cuprate. 

Tc materials. In LaNi02, however, Ni 3d states lie 
just below the Fermi level, with O 2p states located 
3—4 eV below the center of Ni bands. Therefore, p-d 
hybridization, which plays a crucial role in the elec- 
tronic structure and superconductivity of CaCu02, 
becomes much weaker. 



V. SUMMARY 

Aside from the formal similarity to CaCu02, the 
interest in LaNi02 lies in the occurrence of the un- 



usual monovalent Ni ion. As we have found and 
in apparent agreement with experiment, this com- 
pound is a metal, and the "charge state" of a tran- 
sition metal atom in a metal usually has much less 
significance than it is in an insulator. It may be be- 
cause the compound is metallic that it is stable, but 
in this study we are not addressing energetics and 
stability questions. 

Hay ward et a/.0| had already suggested that the 
experimental findings could arise from reduced co- 
valency between the Ni 3d and O 2p orbitals, and 
the 30% smaller value of the hopping amplitude t 
indeed reflects the smaller covalency, as does the 
increased separation between the Ni 3d and O 2p 
bands. It is something of an enigma that in CaCu02 
and other cuprates, LDA calculations fail to give the 
observed antiferromagnetic states, while in LaNi02 
LDA predicts a weak antiferromagnetic state when 
there is no magnetism observed. In the cuprates 
the cause is known and is treated in a reasonable 
way by application of the LDA-I-U method. In 
this nickelate, application of the LDA-I-U method 
does not seem to be warranted (although novel be- 
havior occurs it if it used). Rather, the predic- 
tion of weak magnetism adds this compound to the 
small but growing number of systems (ZrZn2,"25| 
Sc3ln,[2^ and Ni3Ga,_2^J for example) in which the 
tendency toward magnetism is overestimated by the 
local density approximation. It appears that this 
tendency can be corrected by accounting for mag- 
netic fluctuations. 01 
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